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Abstract— This paper studies precipitable water vapor (PWV) in 2015-2016, an El Nino-Southern Oscillation 
episode (ENSO), and causes drought in Indonesia. Using the Global Positioning System (GPS), we obtained the 
PWV value. Expected can detect ENSO phenomenon through PWV value. We processed the GPS observation 
data by Precise Point Positioning (PPP) method using the goGPS software. The obtained values of Zenith 
Troposphere Delay (ZTD) and PWV were validated by comparing with the ZTD values from four IGS stations 
(BAKO, PIMO, NTUS, ALIC), i.e., there was a mean difference of -3.48, the average standard deviation is 13.36 
mm, and the correlation coefficient was 0.94. The variation between the IGS ZTD and PPP ZTD estimates is 
within the measurement error. Our results can be used to study PWV behaviors for meteorological and 
climatological purposes in Indonesia. The values of PWV obtained from station GPS CMAK and CBIT were 
compared with PWV by radiosonde observations. The results showed a good agreement with the coefficient 
correlation is 0.86 at CMAK Station and 0.85 at CBIT station. The PWV values at CBIT stations varied between 
37.08 and 76.01 mm with an average of 57.05 mm. For CTOL, CMAK, and CKEN stations, their average PWV 
is 61.17, 50.84, and 51.75 mm. The standard deviation (STD) value of the PWV at CMAK stations was 11.79 
mm, higher than the other stations. During the El Niño event, the correlation between SSTa and PWV became 
weak at the four stations. 
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I. INTRODUCTION 


The Global Positioning System is the oldest Global 
Navigation Satellite Systems (GNSS). Indonesia's GPS 
network has been used only for positioning, surveying, and 
mappings, such as inland administration, mining, and 
transportation sectors. It has not been utilized to observe 
the Earth's atmosphere, including meteorological and 
climatological studies. In the early 90s, (Bevis et al. 1992) 
first proposed GPS meteorology, i.e., we can use ground 
GPS stations to study water vapor behaviors in the lower 
atmosphere (the troposphere) called precipitable water 
vapor (PWV). 


PWV is an important parameter that has an essential role 
in weather and global climate. One of the factors 
influencing climate change is the ENSO phenomenon, El 
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Nino, Southern Oscillation, and the phenomenon opposite 
ENSO (LaNijfia). We used the four indicators to detect the 
ENSO phenomenon, the sea surface temperature anomaly 
(SSTa), sea surface pressure anomaly (SLPa), and 
multivariate ENSO index (MEI), and rainfall gradient 
(Curtis and Adler 2000). 


This study investigates the behavior of PWV in 2015- 
2016. Based on data from the National Oceanic and 
Atmospheric Administration (NOAA), El Niño in 2015- 
2016 was more robust than El Niño in 1997-1998. During 
El Niño, many parts of Indonesia experience drought. We 
address the question of whether GPS-PWV can detect 
EINifo, especially when El Niño grows and decays. These 
techniques are cost-effective because we use pre-existing 
infrastructure, the GNSS network in Indonesia. In this 
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research, we will focus on ENSO analysis, especially El 
Nino, by GPS PWV monitoring on the island of Sulawesi. 


Il. DATA DAN METHOD 
PWV Retrieved from GPS 


The measurement concept of PWV with GPS is illustrated 
in figure 1. In the model, a microwave signal propagates 
from the GPS satellite to the receiver on the ground, which 
deforms with the seawater's loading mechanisms. GPS 
signals that propagate through the Earth's atmosphere 
suffer from the variability of the ionosphere's refractive 
index and troposphere. Delays are caused by two factors, 
i.e., the excess path caused by the bending of the 
microwave path and the delay along the way caused by the 
refractive index (Suparta 2012). The tropospheric zenith 
total delay (ZTD) consists of zenith hydrostatic delay 
(ZHD) and zenith wet delay (ZWD). ZWD can be 
calculated based on the Hopfield model [H. Hopfield, 
1969]. ZHD can be calculated using the Saastamoinen 
model (Saastamoinen, 1972).We used the Vienna mapping 
functions (VMF1) to convert zenith delays to those along 
the line-of-sight (Suparta et al., 2011).ZWD is calculated 
by subtracting ZHD from ZTD. ZWD was later 
transformed into approximate PWV (in mm) by employing 
the surface temperature measured at a particular site. To 
calculate the PWV, we need meteorological data, such as 
pressure (P), temperature (T), and relative humidity (H). 
Meteorological data are obtained by meteorological 
sensors installed on the GPS Stations to get the accurate 
PWV. 


pets 
baal 


Fig. 1 Propagation of GPS signals to a receiver on the 
ground that covers the sea surface temperature influences 


The conversion from ZWD to PWV is done using the 
formula by (Bevis et al. 1992), 


PWV =x(Tm) ZWD (1) 
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Where 2(Tm) is the conversion factor that varies with the 
local climate and depends on weighted mean temperature 
(Tm, in Kelvin). 


Tm = 0.83663Ts +48.103 (2) 


Where Ts is the surface temperature (in Kelvin), note that 
the Tm equation (2) was obtained empirically by linear 
regression for 15 selected radiosonde stations over the 
Western Pacific (Bevis et al. 1992) 


Dataset 

GPS data for the meteorological study were collected from 
the four stations in Sulawesi Island, Indonesia. We 
selected the El Nifio episode that occurred in 2015/2016. 
We processed the data to study the response of PWV to 
ENSO. The four GPS stations are Makasar (CMAK), 
Kendari (CKEN), Toli-Toli (CTOL), and Bitung (CBIT) 
observatories. Data from the meteorological sensors were 
used for the individual GPS stations. All GPS stations 
recorded the carrier phase with an interval of 30 seconds. 
The stations are shown in Fig. 2. To link the activities of 
ENSO with PWV response, SSTa Oceanic 


Nino Index (ONI) on the Nino 3.4 area was obtained daily 


from the National Oceanic and Atmospheric 


Administration (NOAA). 
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Fig. 2 Location 4 station for study 


Data Processing 


Ground-based GPS data can be processed to determine 
ZTD in two different ways : (1) Baseline network solution 
using double-difference observable (Bevis et al. 1992) and 
(2) Precise Point Positioning (PPP) using indifference 
observable (Byun and Bar-Sever 2009). 


The baseline approach is the most commonly used in GPS 
processing. It is done with two or more receivers 
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simultaneously observing GPS satellites. This technique 
provides precise ZTD estimates, but the separation 
distance between adjacent receivers limits its effectiveness. 
On the other hand, PPP can be performed for a set of 
single GPS receivers to estimate the precise values of 
ZTD. It should be emphasized that in PPP, precise satellite 
orbit and clock correction data are needed. 


In this research, the ZTD results derived using PPP 
processing methods were compared and evaluated with the 
ZTD from the International GNSS Service (IGS) final 
troposphere. One week of GPS observations was collected 
from 4 IGS network stations, 1.e., BAKO, NTUS, PIMO, 
ALIC. This validates that the methods we used to calculate 
the PWV in subsequent studies. For more information on 
the IGS final 
http://igscb.jpl.nasa.gov/. 


troposphere product, refer to 


To process and analyze the data signal and GPS 
meteorology in RINEX format, we used the software 
package goGPS in the PPP model. goGPS is an open- 
source software application developed by (Realini et al. 
2012) since 2007 at Geomatics Laboratory of the 
Politecnico di Milano, Como Campus. Initially, it was 
created in MATLAB but was recently converted to Java to 
expand its users. They started to provide it as a service 
through the web. 


To validate the PWV values from GPS, we compared them 
with those from radiosonde observations around the GPS 
stations we used. Radiosondes provide various quantities 
concerning the atmosphere, including temperature, 
pressure, and humidity. They are often used in validating 
GPS-PWV data as a reliable source of independent data. 
Although radiosondes provide meteorological observations 
with a high vertical resolution, their temporal and spatial 
resolutions are not high. Therefore, we compare GPS- 
PWV values observed at times when the radiosondes are 
released. The PWV from radiosondes is considered to be 
accurate to ~1.2 mm (Feng, Bai, and Williams 2001) 


In this study, we used two radiosonde sites on Sulawesi 
Island operated by the Indonesia Agency of Meteorology, 
Climatology, and Geophysics. They were found to be close 
to two of the GNSS stations used in this study. These 
stations located at the Sultan Hasanudin Airport in 
Makasar (WAAA) and the Sam Ratulangi Airport near 
Bitung (WAAM). The nearest GNSS stations to these 
radiosonde sites are Makasar Observatory (CMAK) and 
Bitung (CBIT). 
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HI. RESULTS AND DISCUSSION 
Troposphere Zenith Delay 


The mean, standard deviation, and the correlation between 
the ZTD estimate obtained from goGPS software with PPP 
processing and ZTD from IGS for the stations over one 
week are listed in Table 1. The mean difference between 
the two processing methods is more petite than -3.48 mm. 
The mean, standard deviation is 13.36 mm, and the mean 
correlation is 0.89. This shows good agreement results in 
ZTD goGPS with ZTD from IGS. 


Table 1 Statistical compilation ZTD processing by goGPS 





software and ZTD from IGS 
Station ID Diff imm] | STO (mm Correlation 
BAKO -3.77 13.58 0.87 
PIMO 6.47 13.60 0.95 
_ ALC | 0.40 | 13.20 pag 
Average | 3,46 | 13,36 ERE 
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Fig. 3 One-week comparison between the IGS final 
tropospheric product with PPP processing solutions at 
BAKO GPS station 


The ZTD estimates for the BAKO station from the IGS 
and our solution are compared in Fig. 3. The results from 
both are consistent, and the deviations between the two 
solutions are well within the level of 1—2 cm. These results 
conform to the findings in (Bevis et al. 1992),(Gao and 
Chen 2004), (Gendt, Reigber, and Dick 2001), and 
(Rocken et al. 1993). An error of 1—2 cm in ZTD equals an 
error of 1-3 mm in PWV, which is considered 
insignificant for this study. Therefore, we use this PPP 
solution in the following analysis, unless stated otherwise. 


Validation of GPS-PWV with 
radiosondes. 


ground-based 


One-month of GPS-PWV estimates at these stations were 
calculated using measurements from the nearest weather 
stations. The GPS-PWV and radiosonde-PWV time series 
at these stations are illustrated in Fig. 4 


The correlation coefficient shows good agreement between 
the two data sets, particularly at the CMAK-WAAA 
station pair. The GPS receiver and radiosonde station are 
collocated. We use these PWV values to study the link 
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between PWV and climate changes, e.g., the El Niño 
episodes. 


PWV and SSTa Variability 


Figure 5 shows the variation of daily averages of the PWV 
from GPS measurements in 2015/2016 at four selected 
CBIT, CTOL, CMAK, and CKEN. This figure shows that 
the PWV at the four stations shows different patterns. At 
the CBIT station, PWV varies between 37.08 and 76.01 
mm, with an average of 57.05 mm. For stations CTOL, 
CMAK, and CKEN, their average PWV is 61.17, 50.84, 
and 51.75 mm, respectively. The standard deviation (STD) 
of the PWV at the CMAK station, 11.79 mm, was higher 
than other stations. 





Fig. 4 Comparison between the GPS-PWV (solid line) and 
radiosonde-P WV for CBIT (upper) and CMAK (lower) 
stations, respectively. The red circles show radiosonde 

estimates of water vapor. 


From Figure 5, PWV at CMAK and CKEN stations show 
a clear pattern that PWV values are higher during the rainy 
season (Oct-Mar) and lower during the dry season (Apr- 
Sept). This tendency was also seen at the CBIT and CTOL 
stations, but there was not so strong. The PWV at CBIT 
and CTOL increase during April-May 2016. Significant 
PWV fluctuation occurred during 2015 at the CMAK and 
CKEN stations in July when ENSO activity disturbed the 
two stations’ atmosphere. 
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Fig. 5 The daily average of PWV variation from GPS 
measurements during 2015 and 2016 


Intense ENSO episodes were recorded in 1982/83, 
1997/1998, and 2015/2016. Such ENSO episodes play an 
essential role in generating extreme conditions in the 
Pacific Ocean and the region around Sulawesi. There are 
two approaches to identify the occurrence of ENSO. In the 
first definition, an ENSO event is considered to occur 
when the 5-month running average of the SST anomaly 
(SSTa) in the Niño 3.4 region is above 0.4 degrees for six 
months (Trenberth 1997). As the second definition, the 
Japan Meteorological Agency (JMA) uses the oceanic of 
Nino 3 (4N-4S and 150-90W). They consider it an ENSO 
event when the 5-month running average of the SST in this 
area (SSTa) exceeds the thresholds of 0.5 degrees or more 
for six consecutive months (Bove et al. 1998). 
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Fig. 6 Scatter plot between GPS PWV and SSTa for Index 
Niño 3.4 region for the year of 2015/2016 


Figure 5 shows the identification of ENSO, according to 
NOAA, based on the first definition. The period of more 
than~12 months shown in red Five to an El Niño episode 
(SSTa> 0,5 ), while the part is shown in blue (SSTa <-0.5), 
is a La Niña event. The gray color part is neutral. From 
this figure, the El Niño event was from the middle of 
January from 2015 to May 2016. Its peak is around 
November 2015, where the SSTa anomaly in the Nifio 3.4 
region 1s 3.28 degrees. 


Monitoring of PWV Variability from GPS During El 
Nino Event 


The correlation analysis was done to show PWV responses 
to El Niño events, i.e., we compare the PWV values 
among the periods of El Niño, neutral, and La Niña. 
During the El Niño episode, its intensity (ENSO index) 
varies in time. We also analyze the PWV changes 
associated with the difference in the index. The responses 
seem to go among the four stations. 


Table 2 The correlation coefficient between GPS PWV and 
SSTa during the El Nino event for the case of 2015/2016 





Station (City) Niño 3.4 Normal LaNina 3.4 





CBIT (Bitung) -0.1606 -0.3460 -0.1469 


CTOL (Toli-toli) 0.1126 0.0748 0.0882 








CMAK (Makasar) -0.1627 0.3395 -0.2384 
CKEN (Kendari) -0.1378 0.4710 -0.0579 
Average -0.0871 0.1348 -0.0887 


Below is the classification of correlation coefficients, 





1. 0. “no correlation.” 

2. 0.20 -0.39 "weak." 
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3. 0.40 -0.59 "moderate" 
4. 0.60 -0.79 "strong." 
5. 0.80 -1.0 "very strong" 


During the El Niño episode shown in Figure .6a., the 
correlation between SSTa and PWV is very weak at four 
stations, 1.e., below 0.20. For neutral conditions (Figure 
6b). The correlation between SSTa and PWV is weak at 
CBIT, very weak at CTOL, soft at CMAK, and moderate 
at CKEN (0.4710). During LaNiña, a weak correlation is 
found only at Makasar station (CMAK). It is very weak at 
the other three stations. 


The weak response to El Niño at these sites may be due to 
seasonal changes from the rainy season to the dry season. 
The relationship between PWV from GPS and SSTa for 
Niño 3.4 did not show any significant correlation, either. 


During the El Niño intensity increase in April-November 
2015 (Figure 7a). Of the four existing stations, a 
correlation between GPS PWV and SSTa is very weak at 
CBIT, weak at CTOL, moderate at CMAK, and soft at 
CKEN. The overall average correlation is very weak. 


tal 
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Fig. 7 Scatterplot between GPS PWV and SSTa from Nino 
3.4region during El Nino event for 2015/2016 (a) the 
decreased intensity and (b) the increasing intensity. 
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Fig. 8 Graph of average PWV value under El Niño 
condition, neutral condition, and LaNiña, at four stations 
in 2015/2016 
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On the other hand, decreasing the ENSO index in 
November 2015—June 2016 (Figure 7b) shows a moderate 
correlation with GPS PWV at CBIT, the weak correlation 
at CTOL, the very weak correlation at CMAK, and CKEN. 
The overall average correction between the declining El 
Niño intensity and PWV is soft. 


Figure 8 shows that the average PWV during LaNifa is 
generally higher than in El Nifio and Neutral conditions at 
all four stations. More statistics details are given in Table 
2. However, further and more in-depth research is needed 
to address the response of PWV to ENSO. 


IV. CONCLUSIONS 


We studied the variability of GPS-PWV and SSTa during 
the 2015/2016 period covering El Niño, neutral. La Niña 
episodes using data from four stations in Sulawesi, 
Indonesia, i.e., CBIT, CTOL, CMAK, and CKEN. The 
GPS-PWV and SSTa showed a negative correlation. The 
decrease in PWV also correlates with the dry season 
lasting from April to September, when water vapor in the 
atmosphere offers a minimum. The selected GPS receivers 
in the Sulawesi region may not be optimal for studying 
PWV changes by ENSO since the correlation between 
PWV and SSTa was weak throughout the period. 
Nevertheless, by comparing the average PWV values in 
the El Niño, neutral, and La Niña conditions, We could 
show that average PWV offers smaller deals in the El Niño 
period than in the La Niña period. Further research is 
needed to find out more details. 
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